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Abstract 

We consider two-nucleon scattering close to threshold. Partial-wave amplitudes are 
obtained by an analytic extrapolation of subthreshold reaction amplitudes calcu- 
lated in a relativistic formulation of chiral perturbation theory. The constraints set 
by unitarity are used in order to stabilize the extrapolation. Neutron-proton phase 
shifts are analyzed up to laboratory energies Ti a b c± 250 MeV based on the next-to- 
next-to- next-to-leading order expression for the subthreshold amplitudes. We find a 
reasonably accurate description of the empirical S- and P-waves and a good conver- 
gence of our approach. These results support the assumption that the subthreshold 
nucleon-nucleon scattering amplitude may be computed perturbatively by means of 
the chiral expansion. The intricate soft scales that govern the low-energy nucleon- 
nucleon scattering are generated dynamically via a controlled analytic continuation. 



1 Introduction 



The last decade has witnessed an impressive progress towards quantitative so- 
lution of the nuclear many-body problem starting from the underlying forces 
between the nucleons. Rapidly increasing available computational resources 
coupled with modern ab-initio few- and many-body methods as well as renor- 
malization techniques to reduce the many-body model space make it nowadays 
possible, to carry out reliable and accurate nuclear structure calculations for 
light and even medium-mass nuclei. One can therefore directly relate the fine 
properties of the nuclear Hamiltonian to the spectra and other properties of 
nuclei without invoking any uncontrollable approximations. It is thus of ut- 
most importance for contemporary nuclear physics and nuclear astrophysics, 
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to develop a detailed, quantitative understanding of low-energy interactions 
between the nucleons based on QCD, the underlying theory of the strong 
force. Remarkable progress along these lines has been achieved in the past two 
decades within the framework of (chiral) effective field theory. In particular, 
accurate nucleon-nucleon (NN) potentials at next-to-next-to-next-to-leading 
order in the chiral expansion have been constructed [TH2] . The corresponding 
developments for the three-nucleon force are underway, see Refs. [3]T4] and 
references therein. 

As a complementary approach, NN scattering can also be addressed from 
the standpoint of the dispersion relations. This method has been formulated 
and extensively explored in the sixties of the last century starting from the 
pioneering work by Goldberger et al. [5], see also Refs. [6fTf8] . Clearly, these 
studies were lacking a systematic theoretical framework to treat low-energy 
pion dynamics which is nowadays available in terms of chiral perturbation 
theory (ChPT). We also mention recent work Refs. PfTU] along these lines, 
where the impact of the left-hand cut emerging from one-pion exchange on 
NN scattering is investigated. For recent applications of dispersion relations 
in combination with chiral perturbation theory to pion- and photon-nucleon 
dynamics the reader is referred to Refs. [TT|12|13lfH] . 




In the present paper, we use these ideas as developed in the recent works 
|llJll2jll3Jll5Jll6j . which we call dispersive effective field theory (DEFT), to 
propose a novel approach to NN dynamics. It relies on the knowledge of the 
analytic structure of the NN S'-matrix and makes use of the chiral expansion 
for the scattering amplitude, which is assumed to be valid in some region below 
threshold, coupled with the conformal mapping techniques to perform reliable 
extrapolations to higher energies. As will be shown below, our approach yields, 
in particular, new insights into the relative importance of contributions associ- 
ated with the two-pion exchange cuts. These results are complementary to the 
ones obtained in earlier studies based on the potentials derived within ChPT 



Our manuscript is organized as follows. In section [21 we briefly describe the 
theoretical method we are using to compute the NN scattering amplitude. 
Chiral expansion for the discontinuities across the left-hand cuts is discussed 
in section [3j The results for NN phase-shifts are presented in section HJ Finally, 
our conclusions are summarized in section [5j 



2 An analytic continuation of the NN scattering amplitude 

A NN partial-wave amplitude has the well-known singularity structure as a 
function of the complex variable s. On the physical sheet, apart from the pole 
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Fig. 1. Cuts of the partial wave NN amplitude in the complex s-plane. The first, 
second and third left-hand cuts start at s = 4m 2 N — (nM^) 2 , n = 1,2,3. The first 
and second right-hand cuts due to inelastic channels start at the one- and two-pion 
production thresholds s = (2itin + nM^) 2 with n = 1,2. Further, hm denotes the 
matching point as explained in the text. 

corresponding to the deuteron bound state in the 3 Si — 3 -Di channel, there are 
right- and left-hand cuts as shown in Fig. [IJ The right-hand cuts correspond to 
intermediate states in the s-channel and start from the two-nucleon threshold, 
s = 4:m%. 

The discontinuity of the amplitude across this cut is determined by the uni- 
tarity condition 

AT(s) = i (T(s + it) - T(s - it)) = T{s + it) p{s) T(s-ie) + ... , (1) 

where the dots stand for inelastic channel contributions which start at the 
pion production threshold, s = (2mjv + M n ) 2 . The phase-space function p(s) 
and the amplitude T(s) in Eq. (0Q) turn into 2x2 matrices for the coupled 
partial waves. 

It is advantageous to discriminate between the two sources of branch cuts in 
the scattering amplitudes. The so-called generalized potential U(s) collects 
all contributions emerging from the left-hand cuts only. Given the generalized 
potential, the full scattering amplitudes may be reconstructed in terms of the 
non-linear integral equation 

t (s)=u{s)+ r ( 2) 

hm\ 7T S' - [M 2 M S' - S - it 

Obviously, the solution of this equation recovers the right-hand cut in agree- 
ment with the elastic unitarity constraints. Since we neglect the contributions 
from inelastic channels in Eq. (JT|), U(s) is well-defined only up to some s = A s 
below which the inelastic contributions are small. We choose this value to 
correspond to the two-pion production threshold A s = (2m^ + 2M n ) 2 . The 
dependence of our results on a particular choice of A s will turn out to be 
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rather weak, see section HI At the matching point s = fi 2 M , the scattering 
amplitude T(s) per construction equals the generalized potential U(s). We 
choose the matching point to be in the middle of the region between the uni- 
tarity s-channel cut and two-pion t-channel cut, i.e. h\l = 4m^ - 2M 2 . The 
proper choice of the matching point is essential for the effectiveness of our 
approach. Our central assumption is a perturbative nature of the scattering 
amplitude near the matching point. Consequently, the generalized potential 
can be calculated at the matching point in the framework of ChPT. Upon 
performing an analytic continuation in the s-variable by employing a suitable 
chosen conformal mapping technique, the generalized potential U(s) can be 
used in Eq. (J2J. 

Consider the left-hand singularity structure of the generalized potential. The 
physically relevant cuts correspond to the t-channel one-pion and several-pion 
exchanges as visualized in Fig. [TJ These cuts start at s — Am 2 N — (nM n ) 2 , 
where n is the number of exchanged pions. We follow here the strategy of [11] 
and split the generalized potential into two parts, 



where Z7inside(s) is calculated explicitly whereas [/outside (s) is analytically con- 
tinued. In Fig. [1] we introduce a domain Q, which defines the separation of 
the 'inside' and 'outside' contributions to the generalized potential. The left 
and right-most point of the domain Q are specified by A 2 (defined below) 
and A 2 . The discontinuity across the left-hand cut from s = 4m 2 N — M 2 to 
s = 4m% — 4M 2 is unambiguously determined by the one-pion exchange con- 
tribution in terms of the nNN coupling constant. It lies inside the domain 
fl as shown in Fig. [TJ The next part of the cut from s = Am 2 N — 4M 2 to 
s = 4m 2 N — 9M 2 receives, in addition, the contributions from the two-pion 
exchange processes. We assume that in this energy region such contributions 
can be computed in ChPT and, therefore, they are also considered as part of 
[/inside (s)- For s < 4:m 2 N — 9M 2 , the two-pion-exchange cut can no longer be 
reliably computed in perturbation theory. This reflects the nature of the tttt 
interaction in this kinematical region with its prominent and non-perturbative 
scalar-isoscalar channel. Therefore, we do not include this part into U\ ns \& e (s) 
and identify 



where AT(s) is the discontinuity of the amplitude across the cut. By definition, 
the 'outside' part of the potential, Z7 utside( s )i has singularities only outside 
the region Q of Fig. [lj It incorporates, in particular, effects of three pion 
exchanges that are technically more challenging to compute. Since U outsi d e (s) 
is analytic inside Q, one may approximate it by a Taylor expansion around 
the NN threshold. The expansion coefficients would resemble the infinite tower 



U(s) = t/inside (s) + [/outside (s), 



(3) 
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of counter terms appearing in the chiral Lagrangian. However, such a series 
would converge only inside a circle with the radius R = 9M 2 and, therefore, the 
generalized potential will be not suitable for using as input in Eq. (J2J) in order 
to compute the partial-wave scattering amplitudes in a controlled manner. We 
recall that the generalized potential is needed at least up to s ~ 4m|r+8mjvM 7r 
which is roughly 6 times farther away from the NN threshold than accessible 
by that Taylor expansion with convergence radius R = 9M 2 . 

Following [TT|23|ll5j . we approximate ^outside (s) in the domain Q of Fig. [I] in a 
systematic way by means of an appropriate conformal transformation. A func- 
tion £(s) having inverse that maps Q onto a unit circle with £(4m^) = is 
constructed. The particular form of the conformal mapping is irrelevant. The 
potential {/outside (s) can now be expanded in a Taylor series in the conformal 
variable, 



KmtsidcO ) = J2U k f (S 
k=0 



k TT d fc [/ ou t s ide (C(£)) 



U, 



k\d£ k 



(5) 



that can be proven to converge in the full domain Q of Fig. [TJ The coefficients 
Uk are determined by derivatives of {/outside at the expansion point s = 4m 2 



n- 



If the expansion in ([3]) is truncated at order n, the first n-derivatives need to 
be computed. We choose the number of terms n + 1 in Eq. (JSJ) to be equal 
to the number of local counter terms contributing to the considered partial 
wave in ChPT at a given order so that there is a one-to-one correspondence 
between Eq. (jSJ) and the chiral expansion. 

To be specific, we use here the conformal mapping suggested in [11] 

= oi (Al-sf-1 = 1 = 1 

(a 1 -2a 2 )(AZ-s) 2 + l , 1 (A 2 - 4m 2 v ) 2 ' 2 (A 2 - A 2 ) 2 ' 

(6) 

A useful property of the transformation (jB]) is that it permits a smooth ex- 
trapolation of the generalized potential to a constant above s = A 2 . Replacing 
in (JSJ the conformal map £(s) by £>(s) defined as 

,f(s) for s < A 2 , , _ 

m 2 ) = o, (7) 

e(A 2 ) for s > A 2 , 

the outside potential is smoothly extended to a constant at energies s > A 2 . 

The final step in obtaining the amplitude in the physical region is to solve the 
integral equation (j2J) with respect to T(s) for a given U(s). A standard way of 
doing that is the so-called NJ D technique [24] . The amplitude is represented 
as 

T(s) = D- 1 (s)N(s) 1 (8) 
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where D(s) has no singularities but the right-hand s-channel unitarity cuts. In 
contrast, the branch points of N(s ) correspond to those of U(s). The unitarity 
condition implies 

J4m 2 N 7T S' — /% S' — S 

The non-linear equation (J2J) reduces to the linear one for N(s), 

at(.) = + £ g "WMFW-uM , (10) 

which can be solved using standard numerical techniques. 

It should be emphasized that apart from the dynamical singularities in the 
complex s-plane discussed above, the partial-wave amplitudes may have kine- 
matic singularities or obey kinematic constraints. In particular, the standard 
JLS amplitudes possess i/i-type singularities and obey constraints at thresh- 
old, i.e. at y/s = ±2rriN. Although, apart from constraints at y/s = 2ttin, they 
are located rather far away from the physical region, it is preferable to remove 
all of them if possible. A general method to achieve this goal for a system 
of two interacting spin- 1/2 particles is presented in [25] and some particular 
cases were discussed in :(ij . We choose a set of amplitudes (which is generally 
speaking not unique) for the NN system as follows. For the uncoupled partial 
waves with J = L and the 3 Po channel, the S'-matrix elements S JLS (s) or, 
equivalently, the phase shifts S J (s) are related with the amplitude T JLS (s) 
which is free of the kinematical singularities via 

S JLS (s) = e 2lSJLS ^ = l + 2z p L (s) T JLS (s) , (11) 

where the phase-space function is defined according to 




For coupled partial waves, one needs, in addition, to perform a linear trans- 
formation of the L = J ± 1 basis leading to the following relation in terms of 
the Stapp parametrization [26] 



( cos2e J e 2i5J " 1 % sin 2e J e i ( (5J " 1+5J+1 ) 
\i sin 2eje^ J - 1+5J+1 ) cos 2eje 2isJ+1 

l + 2ii-fi-^] XFT J U i (13) 
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with the transformation matrix 



U 



2(J+l)m N 




J+l 
J 



,2 



2m jy 

1 

1 ' 



which implies that the NN phase-space distribution has the form 

2 



P 



2 J+l 
8nJ 



1 



A' 



I _|_ 4m N J+l 



J+l ( 1 



V 



J+l / 1 



1 n \ "in 



(14) 



}N_ \ rn N 
s ) Vs 



(15) 



3 Chiral perturbation theory for the left-hand cuts 



We now discuss the derivation of the generalized potential in ChPT up to 
order Q 3 with Q ~ M n referring to a generic soft scale. The following terms 
in the effective Lagrangian are relevant for our calculation: 

C m t = N V (r • (tt x (drf)) N+^-N l5 Y(r- (drf) N 

-^M 2 n N(n-7f)N+^N(d,n)-(d^)N 

-^Na^(r- ((d,n) x (d^))) N . (16) 

Here, tt and N refer to pion and nucleon fields, f denote the isospin Pauli 
matrices while F„ and g& are the pion decay and the nucleon axial vector 
constants, respectively. Further, q are the low-energy constants (LECs) ac- 
companying the subleading pion-nucleon vertices. For more details on the 
effective pion-nucleon Lagrangian the reader is referred to Refs. [2TP5] . 

We apply the effective Lagrangian introduced above to compute the NN scat- 
tering amplitude to order Q 3 . As already pointed out before, our central as- 
sumption is perturbativeness of the scattering amplitude in the region s ~ \j? m . 
We, therefore, use here the standard chiral power counting based entirely on 
the naive dimensional analysis. Although we employ the manifestly Lorentz- 
invariant form of Ci nt without performing the 1 /m^-expansion in order not 
to distort the discontinuity structure across the left-hand cuts, we apply the 
standard power counting rules of the heavy-baryon formulation and keep only 
those diagrams which do not vanish in the heavy-baryon framework at a given 
order. The Feynman diagrams emerging at various orders in the chiral ex- 
pansion are depicted in Fig. Notice that the shorter-range contributions 
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Fig. 2. Feynman diagrams contributing to the NN scattering amplitude at order 
Q % in the chiral expansion. Solid dots (filled squares) denote the lowest-order vertices 
in the first line of Eq. (fTBI) (subleading vertices oc Cj). The ellipses denote diagrams 
leading to shorter-range contributions. 



emerging from graphs which are not shown explicitly in Fig. [2] including loop 
diagrams with contact interactions, 3- and 4-pion exchange diagrams do not 
generate left-hand cuts inside the f2-domain and, therefore, can only influence 
the normalization of the amplitude at the matching point. For the S- and 
P-waves, these effects are absorbed in the corresponding counter terms which 
are available at order Q 3 . For higher partial waves one, in principle, needs to 
evaluate such diagrams explicitly, at least at the NN threshold. We remind 
the reader that at this order, the generalized potential in partial waves with 
L + V > 2, where L (L') denotes the initial (final) orbital angular momen- 
tum, contains only the constant term Uq in the ^-expansion. Given that the 
explicit evaluation of Uq in high partial waves is technically involved and that 
one only expects a small effect of this term in high partial waves due to its 
shorter-range nature, we neglect it in the present analysis. 

Given explicit expressions for the perturbative scattering amplitude, the gen- 
eralized potential can be obtained using the framework outlined in the pre- 
vious section. The explicit calculation of two-pion exchange diagrams in the 
covariant framework is carried out e.g. in Refs. [291130] . For calculating dis- 
continuities it is convenient to employ the Passarino-Veltman reduction of the 
amplitude [3Tf29 32 j in terms of the scalar loop integrals, for which simple dis- 
persion representations are available, see e.g. [29]. Renormalization of the NN 
amplitude requires a special care in the covariant framework. We follow here 
the strategy of Ref. [33J. After applying dimensional regularization in the MS 
scheme, an additional subtraction is done in order to restore the power count- 
ing. The subtracted terms can be absorbed by lower-order counter terms. In 
practice, it is convenient to use the prescription proposed in [29|32] in which 
the subtraction is carried out directly at the level of the scalar loop integrals. 
The details of this calculation will be published separately. 
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4 Results for neutron-proton phase shifts 



We are now in the position to discuss our numerical results for neutron-proton 
phase shifts. In this work, we restrict ourselves to low partial waves where 
non-perturbative effects are most pronounced. In each channel, we solve the 
nonlinear integral equation ([2]) numerically using the N/D method as outlined 
in section [2j The {/i ns id e (s)-part of the generalized potential depends on the 
LECs F^, qa and Cj. Here and in what follows, we adopt the value F n ~ 
92.4 MeV for the pion decay constant. Further, we implicitly account for the 
Goldberger-Treiman discrepancy by using the effective value for g^, 

9a = ^ 1-285 (17) 

m N 

based on g'^ NN /(A'K) — 13.54 [34T35] . Such a replacement is legitimate to the 
order we are working and allows us to take into account the one-pion exchange 
discontinuity exactly. The sensitivity of our fits to the precise value of g n NN 
will be discussed below. By the same reason, we always distinguish between 
the neutral and charged pion masses in the one-pion exchange contribution 
although formally these isospin-breaking effects start at order Q 2 . The correc- 
tions due to the different pion masses in the box diagram appear at order Q 3 . 
We include these corrections perturbatively following the lines of Ref. |36j . 
The impact on phase shifts is, however, almost negligible. Last but not least, 
for the LECs Ci i34 , we adopt the following values from the fit of Ref. [37] : 
Cl = -0.75 GeV" 1 ', c 3 = -4.77 GeV -1 , c 4 = 3.34 GeV" 1 . Interestingly, we ob- 
serve that effects of the order-Q 3 contributions proportional to Cj's are quite 
small for the partial waves considered. In particular, varying these constants 
within the limited range suggested in the literature does not significantly af- 
fect the quality of the fit. These findings are in line with the ones of Ref. |21j . 
Finally, the coefficients £4 entering the truncated expansion of the {/outside (s)- 
part of the generalized potential are determined by fitting the empirical phase 
shifts of the Nijmegen partial wave analysis (PWA) [38] up to = 100 MeV. 
The behavior of phase shifts at higher energies comes out as a prediction. No- 
tice that the number of free parameters in a given partial wave is determined 
by the number of local NN counter terms in the Lagrangian relevant at a given 
order. At even orders Q n (n = 0, 2) there are counter terms contributing to 
the partial waves with L + V = n. 

Our results for uncoupled channels are shown in Fig. [3] In the 1 5'o partial wave, 
there is one free parameter (two free parameters) at orders Q° and Q l (Q 2 
and Q 3 ). It is well known that an accurate description of the 1 Sq phase shift 
requires the inclusion of the subleading contact interaction in order to account 
for a rather large effective range in this channel. Our results agree with these 
expectation. In particular, we observe a good description of the data at orders 
Q 2 and Q 3 , where two fit parameters are available. Interestingly, including 
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Fig. 3. Phase shifts in uncoupled S- and P- waves calculated at orders Q° (dotted 
lines), Q 1 (dash-dotted lines), Q 2 (dashed lines) and Q 3 (solid lines) in comparison 
with the Nijmegen [38] (filled circles) and SAID [39J (filled triangles) PWAs. 

only the first left-hand cut Am 2 N — AM 2 < s < 4m 2 N — M 2 explicitly in U- msidc 
leads to fits of a comparable quality. On the other hand, the fit appears to be 
rather sensitive to the value of g n NN coupling constant. Treating it as a free 
parameter yields a value for g 2 NN which lies within 10 percent of the empirical 
one if one fits the data up to Ti a b = 100 MeV. This difference gets even smaller 
if the fit is restricted to a narrower energy region around the threshold. A 
similar sensitivity to the value of g n NN is observed in the uncoupled P- waves 
of Fig. El 

For P-waves, there are no free parameters (one free parameter) at orders Q° 
and Q 1 (Q 2 and Q 3 ). The calculated phase shifts at orders Q° and Q 1 (Q 1 ) 
in the 1 P 1 ( 3 Pi) partial wave show a steep rise already at rather low energies. 
This reflects the appearance of nonphysical resonances generated by solving 
the N/D equation. These artifacts disappear at orders Q 2 and Q 3 . In the 3 Pq 
channel we observe a large contribution of the lowest-order counter term even 
at low energies which agrees qualitatively with the findings of Refs. [lOPT] . 
Given that there is just one free parameter in each of the P-wave, the phase 
shifts are remarkably well reproduced at order Q 3 . 

We now turn to the coupled channels, see Fig. HJ and first address the conver- 
gence pattern in the 3 Si- 3 Di partial waves. The lowest-order (i.e. Q°) results 
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show a large disagreement with the empirical phase shifts. In particular, the 
mixing angle e\ even comes with the wrong sign. The order-Q 1 corrections 
emerging from the box diagram, see Fig. [21 are large and lead to a strongly 
improved description of the data. Notice that there is just one free parameter 
at both orders Q° and Q 1 . This suggests that the long-range physics associated 
with exchange of pions plays a prominent role in this channel which agrees well 
with the results based on chiral and phenomenological NN potentials. Taking 
into account the Q 2 - and Q 3 -contributions leads to further improvement (in 
part, on the cost of having two more adjustable parameters). The situation 
in the 3 P2- 3 -p2 is qualitatively similar: one observes a large disagreement with 
the data at order Q° which is strongly reduced by taking into account the cor- 
rections at order Q 1 . Notice that the results at this order are parameter-free 
predictions. The higher-order corrections bring our results into a better agree- 
ment with the data (on the cost of adjusting a single free parameter to the 
3 P2 phase shift). The agreement with the empirical phase shifts at order Q 3 is 
reasonably good, especially given that there is only one adjustable parameter 
at this order. Note that because of the channel coupling the fit is performed 
simultaneously for the three phases that sometimes leads to situations when 
although the overall fit improves the results get worse for some particular 
phases when going to higher chiral order. 
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Lab. Energy [MeV] Lab. Energy [MeV] 

Fig. 5. Sensitivity of the 1 Sq and 3 Si- 3 -Di partial waves to the variation of A s 
at order Q 3 . Light shaded bands correspond to the variation of A s in the range 
2m n + M n < A s < 2m n + 2M n , whereas dark shaded bands correspond to the 
range 2mj^ + 2M n < A s < oo. 

Apart from the counter terms, there are several parameters in our scheme 
whose particular values were set by physical arguments. First, the matching 
scale fiM should, of course, be taken far away from the non-perturbative regions 
in the s— and t-channels. That is why we choose /ia/ to lie between s- and t- 
channel two-body thresholds. At sufficiently high order in the chiral expansion, 
the results should become independent on the particular choice of \iu since 
the dependence on \i M can be compensated by local NN counter terms. At 
order Q 3 , we observe a rather week dependence on [i M in most of the partial 
waves we have studied. For example, varying the position of hm within the 
region Am 2 N — 3.5M 2 < [i 2 M < Am 2 N — 1.1M 2 affects the 3 D X phase shifts at 
-Ei a b = 250 MeV, where the effect is the largest, by less then 10%. The second 
source of uncertainty is associated with the choice of conformal mapping and, 
in particular, of the value of A s . To the order we are working, the variation of 
A s only affects the S-waves, where the expansion of £/ outsi d c goes beyond the 
zeroth order in £. Small variations of A s do not lead to sizable effects on the 
results. For an illustration, we show in Fig. the effect of the variation of A s 
within the extremely large range 2m n + M n < A s < oo. Both limits are highly 
unphysical: the lower one is located very close to the energy region we are 
interested in and thus makes the ^-expansion there unreliable while the upper 
one assumes the irrelevance of inelastic cuts in the whole energy region. One 
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can see that even within these extreme limits, the variation of the phase shifts 
is rather modest. Only in the 3 Si phase shift and in the mixing angle e\ one 
observes a sizable variation at higher energies. Varying A s within a reasonable 
range does not affect our results in a significant way. The dependence of the 
results on the parameter of conformal mapping A is much weaker than the 
dependence on A s and will be given in a separate publication. 



5 Summary 

In this work, we performed an analytic continuation of the subthreshold NN 
amplitude into the physical region assuming the validity of ChPT for the 
amplitude at some matching point below threshold. The amplitudes at the 
matching point and the discontinuities across the nearby left-hand cuts are 
calculated to order Q 3 using a manifestly covariant version of ChPT. The 
contributions from distant left-hand cuts are extrapolated to higher energies 
using conformal mapping techniques. The resulting generalized potential is 
used to reconstruct the full scattering amplitude by means of partial-wave 
dispersion relations supplemented by the unitarity constraint. Within our ap- 
proach, the short-range part of the generalized potential is parametrized in 
a systematic way in terms of the Taylor expansion of f/ ou tsidc(s) i n a confor- 
mal variable £(s). From the conceptual point of view, this is analogous to the 
representation of the short-range part of the nuclear force in terms of contact 
interactions with an increasing number of derivatives. We have determined 
the corresponding coefficients entering the expansion of U ou t S ide{s) in S- and 
P- waves by fitting them to the Nijmegen PWA up to the energy of T/ a b = 100 
MeV. The obtained fits and predictions for phase shifts at higher energies are 
in a reasonable agreement with the empirical PWAs up to = 250 MeV. 
The quality of the fit at order Q 3 is comparable with the one obtained in 
conventional approaches based on the chiral expansion of the NN potential. 
In all partial waves, the expansion for the phase shifts seems to converge when 
going from the order Q° to Q 3 . This supports our assumption of the validity of 
a perturbative expansion of the NN amplitude at the matching point. We also 
studied the dependence of the results on the particular choice of the matching 
point hm and the conformal mapping parameter A s . We found our results 
to be rather weakly dependent on /im and A s provided these parameters are 
varied within a physically acceptable range. 

Our detailed analysis provides new insights into the relative importance of 
contributions associated with the one- and two-pion exchange cuts. Such con- 
clusions are more difficult to draw in calculations based on a potential scheme, 
where the left-hand cuts in the scattering amplitude emerge not only from 
the potential but also from its iterations within the dynamical equation. As 
an important outcome of our investigation, we found unambiguously that NN 
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scattering in the 3 Si — 3 -Di channel shows a strong evidence of the lowest-order 
two-pion exchange which in the conventional, Schrddinger-equation-based ap- 
proaches corresponds to the first iteration of the one-pion exchange potential. 
In particular, taking into account the leading two-pion exchange is absolutely 
necessary to achieve a good description of the mixing angle e±. 
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